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Motor Speed Est [mation for stabilized n Iotor contool 



CROSS-REFERENCE TO 
This application is a 
5 fded September 4, 2002 



iJealted APPLJCATIOM 
continutlion-in-partof U.S. AppHcatioh 



FIELD OF THE IN VENTIC|N 
The present invention is 



rela tes to the field of motor speed < oatrol 



I this 



10 BACKGROUND 

Mow controllers govern 
torque, motor shaft position 
operation. Alternating curr< nt 
encoUerless feedback for cobtrol 

15 estimators arc based on the 
measurement of the current |or 
rotation induces other 
construction (rotor bars or 
be generated by detecting 

20 ihem to known mathematical 
speed. Such a known 
Induction Machines", P. 



asf ect 



£ nd 



wt ve; 



SUMMARY 

2S In accordance with one 

estimating a speed of a motor 
wave sensed at the motor 
representing the current 
corresponding number of f irst 

30 the orthogonal pairs being 3rth< 
frequency pairs as ^ orthpgonal 
identifying the frequency 
between the current wave 



operation of a motor in order 
and other such characteristics 
(AC) variable speed (motor 
of the motor. Encodcrles' 
>rtnciple that the motor speed 
voltage waveform at a moto - 
in the currem waveform 
inding gaps) of the motor's ro 

induced speed related 
models reiatmg the induced 

model may be, for 
»r; Gordon and Breach, New 



cani 



frequ sncies 



t icse 



math imatical 



Algei 



the« 



of the present invention 
comprising; determining a 
frequency pairs from a set 
; fittmg components of a 
orchogonal pairs in a set of 
ogonal to the frequency 
pair subsequent to the 
] lairs that provide a reduction in 



pa rs; 



md the set of weighted orthogc nal 



Serial No. 10/241^90, 



0 direct activation, speed, 
(jfthe motor during 

drives (VSD) may have 
speed controls or 

be estimated from a 
connection. Motor 
due to the physical 
or. A speed estimate can 
freq^ lencies and comparing 

1 requencies to the motor 
le/The Nature of 

fork, 1965. 



exc mple 



is provided a method of 
dorrelaiion between a current 
weighted frequency pairs 
motor control signal to a 
kveighted orthogonal pairs, 
; fitting each of the 
1^$t orthogonal pairs; 
the mean squared error 
pairs that satisfies a 
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10 



15 



20 



25 



30 



criteria; and comparing d sn td 
motor speed harmonic mode 



In accordance with another 
of estimating a speed of a 
comprising: determining a 
frequency pairs from a set 
fitting components of a motcfr 
orthogonal pairs in a set of 
orthogonal to the Aeqaency 
pair subsequent to the first 
frequency pairs in a mean 
wei^ted orthogonal pairs 
wave and the first orthogonal 
and the coirespondinj 
deteimined corresponding 
from the determined 
model to determine an 



V eic 



Dg frequ sncy 



1 concsj on( 



[ estm ate 



In accordance with a forthei 
method of estimating a 
current wave sensed at the 
frequency pairs representmj ; 
signal TO a coiresponding m mber 
orthogonal pairs, the orthog[>nal 
comparing a subhannonic 
identify regions in which 
a frequency less than a 
frequencies in the regions 
smallest hamionic of the 
harmonics pair of the corre 
with the harmonics speed 



imot>r 



frequencies from the identifi sd frequency pairs witih a 
to determine an estimation of i he speed. 



repn;seniing1 



a spect of the present invention 
m( tor using a current wave scnsefi 
c< Ttelaiion between a current wal/e 
01^ weighted frequency pairs 

control signal to a corresponding 
ighted orthogonal pairs, the o 
}airs; fitting each of the freque jicy 
ofthogonal pairs; determining a 
error between the curren 
according to a mean squared 
pairs; determining die 
pairs; removing imwantec 
frequency psurs; and comparing 
ding fi^equency pairs with a 
for the speed. 



lemr 



! redttcti ^ns 



aspect of the present invention 



mot }T speed comprising: deteimini ig a correlation between a 
I loior and frequency pairs from a set of weighted 
the current wave; fitting com{ onents of a motor control 
of first orthogonal pairs n a set of weighted 
pairs being orthogonal to i he frequency pairs; 
the current wave with a ha monies speed model to 
a corresponding haimoi ic, the subhaimonic having 
control signal; identifying a laimonics pair of 
ijaving a separation from each o her of no greater than a 
control signal^ wherein on( ( 



f t)ml 



tc locates 



m>tor< 



t^erc is provided a method 
at the motor, the method 
sensed at the motor and 
the cunent wave; 
number of first 
thogonal pairs being 

pairs as the orthogonal 
] eduction for each of the 
wave and the set of 
between the current 

that satisfy a criteria 
frequencies from the 
he remaining frequencies 
I notor speed hannonic 



there Is provided a 



;ponding harmonic; and compa ring the harmonics pair 



11 lodel to determine an estimate 



of the frequencies is the 



or the speed. 



2 
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be ween i 



mp ^nents < 



In accordance with yet 
system for estimating a mote 
determining a correlation 
pairs from a set 

5 noechanism for fitting co] 
number of first orthogonal 
of the frequency pairs as the 
the orthogonal pairs being 
mechanism for detcrminingja 

10 wave and the set of wei 
identifying the frequency 
estimation mechanism 
frequency pairs with a moU|r 
speed; and a controller in 

15 mechanism, the mse reduction 
coordinating tihe process 



L inventii »n 



anotller aspect of the present i 
r speed comprising: a 

a cuiicnt wave sensed at 



of weighted frequency pairs representing th ; 



of a motor control signal 
l|ftirs in a set of weighted orthogpnal 

orthogonal pair subsequent to 
olihogonal to the frequency pair 
reduction in the mean square I 
ighte^l orthogonal pairs; a pair comp Mison 

having a reductim that satii fies 
for cfimparing desired fiequencics f 
speed harmonic model to 
cjkmmunication with the conelat 

mechanism and the speed 
oflestimaiiiig the motor speed. 



I aspc :t 



Lcumnt 
r pail s 



mc tor < 



cun^nt 



In accordance with an 
estimating a motor speed 

20 correlation between a 
of weighted firequency 
Ottiing components of a 
orthogonal pairs in a set 
orthogonal to the fiequenc f 

75 subharmonic from the 
regions in which to locate 
frequency less than a mot< 
for identifying a harmonic^ 
from each other no 

30 wherein one f the frequeni 
speed estimation mechanl 
harmonics pairs with a 
speed; and a controller in 



» gieatei than 



I correlatii »n 



there is provided a 
mechanism for 
the motor and frequency 
current wave; a fitting 
to at corresponding 

pairs and fitting each 
he first orthogonal paurs, 
; an mse reduction 
error between the current 
mechanism for 
a criteria a speed 
"Om the identified 

an estimation of the 
on mechanism, the fitting 
eltimatton mechanism for 



detiaminei 



of the present invention then 
c )mprisjng: a correlation mechai lism 
wave sensed at the motor anil 
representing the current wave ; 
control signal to a correspoi iding 
0^ weighted orthogonal pairs, the 
pairs; a region determination 

wave with a harmonics spefd 
i corresponding harmonic, the 
control signal; a correspondiifg 
pair of frequencies in the two 
a smallest harmonic of th s 
ies m the harmonics pair is the 
;m for comparing desired frequ 
halmonics speed model to deierm 
:ommunication with the correlation 



is provided a system for 

for determining a 
frequency pairs from a set 
a fitting mechanism for 

number of first 
orthogonal pairs being 
nechanism for comparing a 

model to identify two 
iubhannonic having a 

frequencies mechanism 
regions having a separation 
motor control signal, 
corresponding harmonic; a 
ncies from the identified 
ne an estimation of the 
mechanism, the fitting 



3 
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mechanism, the mse reductio i 
cooixiidating the process of e 



mechanism and the speed < 
timating the motor speed. 



In accordance with an aspect 
5 readable medium having stoipd 
estimating a speed of a moto 
wave sensed at the motor an< i 
representing the current wavt 
corresponding number of firjt 
10 tlie orthogonal pairs being 

frequency pairs as the oithoj onal 
identifying the frequency ths 1 1 
the current wave and the set of 
comparing desired firequenci as 
IS harmonic model to deteimir s 



of the present invention there 

thereon computer-cxecutabl 
comprising: determining a 
frequency pairs from a set of 
; fitting components of a moto r 
orthogonal pairs in a set of 
oi^ogonal to the frequency pain ; 
pair subsequent to the firs : 
provide a reduction in the 
wei^ted orthogonal pairs 
fix>m the identified frequency 
an esrimation of the speed. 



aspec t of the present i 



stcredi 



^anl 



In accordance with an 

readable medium having 

estimating a motor speed 
20 wave sensed a the motor 

representing the current wate; 

corresponding number of fi^st 

the orthogonal pairs being 

subharmonic from the currcfit 
25 in which to locate a corresppnding 

less than a motor connt)l s\\ 

a separation from each othcjr 

signal, wherein one of die 

comparing the pair of frequlencies with 
30 an estimate for the speed 



In accordance with an aspc 
readable medium having 



est mation mechanism for 



provided a computer 
I instructions for 
coi relation between a current 
' vcighted frequency pairs 
control signal to a 
ighted orthogonal pairs, 
fitting each of the 
orthogonal pairs; 
squared error between 
satisfies a criteria; and 
pairs with a motor speed 



W51C 



mem 



ttat! 



\ correU tion 



I mot >r 



invention there 
ihereon computer-executabjle 
cofnpnsing: determining a 
frequency pairs from a set of 
; fitting components of a 
orthogonal pairs in a set of 
c(rtbogonal to the frequency pairk, 
wave with a harmonics spee I 
harmonic, the subhaimqnic 
1 j[nal; identifying a pair of frequim 

no greater than a smallest haivionic 
f|equencies in the pair is the 

the motor speed hafmonic 



:t of the present invention there 
st|)red thereon computer-executable 



s provided a computer 
instructions for 
between a current 
veighted frequency pairs 

control signal to a 
eighted orthogonal pairs, 
; comparing a 
model to identify regions 

having a frequency 
icies in the regions havmg 
of the motor control 
spending harmonic; and 
lie model to determine 



is provided a computer 
instructions for 
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10 



15 



20 



25 



30 



orthogonal to the frequency 
pair subsequent to the first o 



estimating a speed of a moto ' using a cunent wave sensed a : the motor, the method 
comprising: determining a ci irrelation between a current wa /e sensed at the motor and 
fiequency pairs from a set oi weighted frequency pairs repn senting the cunent wave; 
fitting components of a moK r contfo] signal to a corresponc ing number of first 
orthogonal pairs in a set of v eighicd orthogonal pairs, the o ' 

)airs; fitting each of the frequci icy pairs as the orthogonal 
thogonal pairs; determining a i eduction for each of the 
frequency pairs in a mean sq iared error between the curreni 
weighted orthogonal pairs ac cording to a mean squared errc r between the current 
wave and the first ortbogonafl pairs; determining the reducti ms that satisfy a criteria 
and the coircsponding frequ mcy pairs; removmg unwanted I 
determined coiresponding fi ;quency pairs; and comparing i ie remaining frequencies 
from the detennined corresp )nding frequency pairs with a ijiotor speed haimonic 
model to determine an estim ite for the speed. 



BRIEF DESCRIPTION OF 
The present invention will b 



the drawings in which: 



FHE DRAWINGS 
; described in conjunction withlt 
Fig. 1 illustrates a motor cor u:oI system according to an em Todiment of the present 
invention; 

Figs. 2A-B illustrate a method of estimating the speed of a ^otor according to a first 
embodiment of die present i ivention; 

Fig. 3 illustrates a speed esti nator of the motor control system according to an 
embodiment of die present i ivention; 
Fig. 4 illustrates Ae frequen :y spectrum in ^ch die speedl estimator detects speed 
related information; 

Figs. S A-B illustrate a methU of estimating the speed of a jnotor according to a 
second embodiment of the p nesent invoition; and 

Fig. 6 illustrates a speed esti nator of the motor control system according to the 
second embodiment. 

DETAILED DESCRIPTIOT I 

Figure I illustrates a motor < ontrol system 100 according td an embodiment of the 
present invention. A vector controller 1 02 in the motor cor crol system 1 00 drives an 



5 
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AC induction motor 104 via 
motor 104 is driven by a 
produced l^y the vector 
from a rotor flux reference 
5 rotor of the motor 1 04 from 
vector controller 102 provid4s 
flow of power from a power 
receives the cono-ol signals 
source 1 1 6 to the motor 

10 

The AC induction motor 10^ 
response to the power provided 
result of a magnetic field 
frequencies representative 
15 during rotation due to the pljysical 
gaps, etc.). This ctirrent is 
1 12 and forms the cunrent vsfive 
estimate of the motor speed 



m inverter 1 10 and motor drive 
knokvn frequency, th fundamental 
connpUer 102. The vector controllc 
, a speed reference 108 and 
speed estimator 114- Based 
control signals 124toti)e 
source 11 6 to the motor 104. 
24 from the vector controller 
104 may be regulated thereby. 



thu 



includes a stator and Oie rotor 
to the motor 104. Current 
pn|duc6d by power supplied to 
the rate of rotation of the rotoi 
construction of the rotor 
ijieasured by sensois 122 at ihe 
from which the speed 



20 The rotor flux reference 1 0( 
reference irame for 
can be used to control the 
current. 



tn 1 



u iesl 



25 The speed reference 108 is 
The vector controller 102 
motor speed obtained from 
signals 1 24 produced by th< 
104 to align the motor speep 

30 frequency, also called die 
produced by the vector 
104 so^at it is the same 



01 



mv *rter 



>5 



1)2- 



Cnot shovm) that rotates in 
iftduced in the rotor as a 
stator includes 
and frequencies produced 
(e.g. rotor bars or winding 
1 Qotor drive connections 
114 derives an 



estin ator 



is a measure of the rotor flux 
measunjment of torque in the motor 
tclrque by regulating the rotor 



ope ating 



indication of the desired 
the speed reference 108 as i 
he speed estimator 114 can 
vector controller 102 are 
with the speed reference 108 
(laiTier drive frequency (CDF), 
controller 1 02 is generated to adjus t 
speed reference 108 



as the I 



connections 1 12. The 
frequency, which is 
102 receives information 
estimated speed of the 
this information the 

110 for control of the 
the inverter 110 
, power from the power 



1 inkage and provides a 
loji. This reference frame 
fluk-linkagc and stator 



(r 



speed of the rotor, 
baseline against which the 
. The control 
to control the motor 
That is, the fiindameTTtal 
the amplitude of the CDF 
the speed of the motor 



be :ompafed. 
gencated 



it 61 32370045 



T-794 P. 020/065 F-051 



A voltage feedback may be ] rovided fiom the motor drive c omiections 11 2 to the 
vector controller 102. 



The vector controller 102 
5 vector controller 1 02 can 
methods, such as the one set 
Oxford Science Publications 
controller similar to the 



provides the speed estimator ] 14 
be )ased on any of a number 
forth in "Vector Control of AC 
1990. For a synchronous 
vector controller 102 may be used. 



10 The speed estimator 1 14 det 
of the current wave sensed 
the CDF 118, Forlowspee( 
receives a selected carrier 



amines an estimate of the mot^r 
qy the sensors 122 at the motor 
operations, the speed estimatcbr 
frequency (SCDF) 1 20 froi i 



diive 



15 The CDF 118 provided to 

produced by the vector controller 
obtained from the motor dri^ 
during current wave obtalniienty 
controller 102 that the curre 

20 In response, the vector contioller 
predetermined amount of tirhe 
wave samples from the 



^senfDTS 



When The CDF values are 
25 fraction of the base speed 
signal 124. The SCDF 
CDF 118. The CDF 118 
related harmonics in the 
speed estimator 1 14, The 
30 have a sufficient 

the motor's efficiency due 
estimator 114 to estimate 
CDF 11 8 are closely space< 



viththeCDF 118. The 
of kijown vector control 
Ndachines " P. Vas; 
moit)r a motor speed 



speed based on samples 
I Irive connections 1 1 2 and 
1 14 alternatively 
the vector controller 1 02. 



carier 



may 



speed estimator 1 1 4 is the 
104 ai the time when the 
^e connections 112. In order to 

the speed estimator 1 1 4 
It at the motor drive connectioiiis 
102 may produce a constant 
during which the speed estii^ator 
122. 



SI aall 4 



; opera ting 



tteni 



and the motor 104 is 
the CDF 118 and the SCDF 
120ls of a higher frequency but 
anb the SCDF 120 are added tog^dii 
cu^ent wave fiom the sensors 
120 is selected such that 
sepamtionjfor detection. Although addin] 
waste of some power, the 
motor speeds where the freq 
together (e.g. I Hz) making it 



:lov rer 



122 to 



tie 



iSOF 



drive frequency 
ruirent samples where 
provide a constant CDF 
inform the vector 
1 12 is being sampled. 
CDF for a 

114 obtains current 



at only a small 
120 may form the control 
amplitude than the 
ler to produce speed 
be detected by the 
speed related hannonics 
the SCDF 120 decreases 
120 enables the speed 
lencies induced by the 
lard to distinguish them. 
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10 



15 



20 



23 



30 



The SCDF 120 may be a prdset 
60 Hz for all occasions thai ihe 



The speed estimator 1 14 use i 
in a Real Time Fast Orthogopal 
representative of the motor 



frequency used regardless Q|rthe speed (e.g. 30 Hz or 
SCDF 120 is used). 



ibe ! 



theCDFll8/SCDF120and 
Search (RTFOS) process 
sbeed in order to produce the m4>tor 



sampled current wave 
tojestimate harmonics 
speed estimate. 



The inveiter 1 10 receives thi ; control signals 124 from the 
control flow of power from fie power source 116 to tfie 
received from the speed estii nator 1 1 4 does not meet certaii 
vector controller 1 02 enhanc es the CDF 1 1 8 by adding the 
combined CDF 1 18 and SCIpF 120 fonn the control signal 
inverter 110- The inverter 1 



that couple power from the | ower source 1 1 6 to the motor ] 



Pigs. 2A-B illustrate a meth< 

1 14 according to an embodii lent of the present invention, 
on the Fast Orthogonal Sean h (FOS) method. FOS uses a 
functions (e.g. a portion of a 



sampled signal (e.g. the curr ;nx wave). FOS searches for 
longer than a specific sampli ng window as well as those 



periods in the sampling win< ow, 



mo or 



When the motor 104 is in 
1 12 from the rotor of the 
102 may add the SCDF to 
The sensors 122 connected 
this current wave. Thevectdr 
wave from the motor drive c 
controller 102 can produce a 



op ^ration < 



The measured current wave 
from the vector controller 1 



vpctor controller 102 to 
motor 104, When the speed 
minimum criteria, the 
$CDF 120 thereto. The 
24 that is provided to the 
0 may be formed of an array o f field effect transistors 

04. 



!Ct< 



4 200 of estimating motor spec d by the speed estimator 

he method 200 is based 
of non-orthogonal 
sine and a portion of a cosine \ ^ve) and fits them to a 

frequencies with a period 
a fractional nimsber of 



Willi 



t) 



current is induced in th( 
104. Ifthe CDF is small 
CDF. thus the SCDF may alsj> 
die motor drive connections 1 
controller 102 is inforaied in 
)nnections 1 12 is going to be 
constant CDF. 



s received from the sensors 1 2! 
op in step 204, The measured 



motor drive connections 
A the vector controller 
be received in step 204. 
1 2 detect and measure 
step 202 that the current 
n easured so that the vector 



and the CDF is received 
ciMrent wave is convened 
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to a digital sampled current 
wave provide a sampled 
performed. A predetermined 
in step 208. Each of the N 

s by a fixed amount of time 
collecting the samples may 
limited by the ability of the 
estimated speed or measure( 
speed regulator) improves 

1 0 second that are measured, 
speed estimator 1 14 then 
rated motor speed will 



15 



20 



vave in Step 206. Ifthesenson 
curfent wave that is already digital 
number of samples N of the 
s[unples may be taken such that 
consecutive samples. The 
( epend on the number of speed 
peed estimator 1 14 to estimate 
samples. Accuracy of the 
an increase in the number o 
10 10 speed estimates per 
control of the motor 104 operatinj 
be stable. 



^ith j 
I/6i 



genei ally 



The currem wave samples cin be modeled as a weighted sub of candidate frequency 
pairs by: 



where Pj^in) = sm 



land 



[gt ted ( 



where y(n) is the sampled ci 
candidate functions p^(n) 
expansion (sum of the weij 
P2in(n) and P2m+i(n) reprcsen^ng 
frequency Fs and the set 
are chosen to represent fieqi 



w y 



25 The sampled current wave i< 
210. The candidate &equcn( 
pairs (e.g. the candidate func|tions) 
remain die same every time 
can be represented by a wei{ 

30 The candidate fipcquency 
sampled current wave can \h 



122 detecting the current 
then step 206 need not be 
c jjment wave are measured 
^ch sample is separated 
i mount of time spent 
stimates per second as 
the speed and store the 

controller 102 (i.e. 
speed estimates per 

are produced by the 
below 2% of die base 



vec or 



sec ond < 



Pj„^,(n) = cos 



I ai d 



I p. 

rrent wave which is a weightcc^ (aj 
the CiTOr ^(n) between y(n) 
candidate functions). The 
the candidate function are 
of candidate frequencies fin. The 
dencies between QHz and the CJ|)F 



correlated to a set of candidau 
pau:s aie chosen to repre 

having particular firequei^cies 
he method 200 is executed 
hted combination of candidate iTunctions 
represent the candidate fiinci 
derived. A subset of candidal ; 



(1) 
(2) 

combination of the 
the functiooal 
andjdate function pairs 
>ased on the sampling 
c< ndidate frequency pairs 



frequency pairs in step 
and cosine function 
of interest and 
sampted current wave 
and an error, 
on fiom which the 
frequency pairs that is 
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20 



25 



noncontiguous can be used 
perfarmed by the speed 



IS 



estu MOT 



Correlation of the current w^ve 
5 perfonned via known correlation 
Fourier Transform (FFT) 
candidate frequency pairs 
filled FFT. Zeros are addedho 
addfaig 1 024 zeros to the en< 
10 Zcro-fiUing increases the 
by the FFT. If his known 
and addition operations in 



with the set of candidate frequency 
techniques or by using zei o-filled, 
the case of the later correlatio 
chosen to correspond to frequ ency 
the sampled signal to mcrea$e 
of 1024 sample values to 
nu[iiber of candidate fiequeiicies 
that the FFT is zero-filled, 
FFT algorithm can be remov^. 



tiiatl 



I tie 



The candidate frequency pal^ are detenniaed to correspond with the FFT bins and to 
1 5 have a higher frequency res< iution than the FFT. The cand ( 

separated from each other b: approxunately four times the ; TT bm size. 



Correlation of the sampled 
be perfonned using a modifj&d 



.Ink 



where k is the frequency bid index, 
samples in the time series ai er 
be expressed as sinusoids anb 



FFT 



The sampled current wave 
frequency bins calculated. { 
bin frequencies and the reso ution 



M 



30 where 0<c<M/2. 



a means of reducing die numt)er of calculations 
114. 



) inert ase 



pairs may be 
tiieci, adjusted Fast 
1 calculation, the set of 
bins of the zero- 
their length, such as by 
the set length to 2048. 
4alculated for correlation 
s{me of the multiplication 



date frequency pairs are 



qurrent wave and the set of canc$daie 
Fast Fourier Transform (FFT) 



n is the time index and 
it has been zero-filled. This 
cosinusolds resuhing in 
(2^ 



A is the number of 
complex equation (3) can 



is zero-filled to M points to 
^ince the candidate frequencies 
oftheFFTisRpn=F,/M 



frequency pairs may 
given by : 

(3) 



(4) 



mcrease the number of 
should coincide with FFT 
hen 

(5) 



10 



04-F9b-13 09:46an) From-KIRBY EADES GALE BAKER 



tl6t32970045 



1-794 P. 024/065 F-051 



10 



IS 



20 



25 



30 



A suitable type of FFT is th< 
The DIF FFT processes *e 
three-quarters of the input is 
by these zero elements need 
computations perfonned by 



While the sampled current 
candidate functions pmCn), ti 
combination of weighted 



real valued FFT with Decimatjon 
nput in the order collected, 
zero, the multiplication and adtiition 
not be performed, thus reducing 
he method 200. 



in Frequ^cy 05IF). 
the bottom half or 
operations required 
I the number of 



V ave 



on logonal 1 



where gm are the weights foi the 
When the weights gm relatin i 
wave y(n) are determined th 
non-orthogonal candidate jftjjictions 
derived from the candidate 
similar to the candidate 



Such derivation of the orthoi 
the Gram-Schmidt ortfaogon&lization 
described in **Fast Orthogon il 
Estimation;* K. Adeney and M. 
Processing. Vol. 141, No. 1 



y(n) can be represented as 
e sampled current wave y(n) 

functions W(n(n) and s^me 



\ an 



a weighted combination of 
also be presented by a 
error £(n} as follows: 

(6) 



orthogonal functions, 
the orthogonal functions Wm(ik) to the sampled current 
y can be used to detcimine the weights atn relating the 
pni(n) to y(n). The orth igonal functions ate 
llmctions and may be represenix d by orthogonal pairs 
freq lency pairs. 



;ofial functions may be perfonried 
algorithm. The Gram 
Search for Array Processing 
Korenberg; FEE. Proc. Vlsijal 
•eb. 1994, herein incorporated 



using a modification of 
Schmidt algorithm is 
md Spectrum 

Image Signal 
jy reference. 



If the CDF is not less than 1 1z or the motor speed is greats 
as determined in step 212 thi in the CDF is fit as tiie first can dldaze 



the orthogonal pairs set in st ^p 214, 



The process of fitting the ca^didace 
orthogonal to all previously 
As the orthogonal fbnciions 
detennuied as part of the fitting 



frequency pairs involve 
itted functions based on the ca ididate 
derived, the weights of thcs 
process. 



than 2% base rated speed 
frequency pair in 



deriving a function 

frequency pair, 
orthogonal iunctjons are 



11 
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An orthogonal function can 
minus the wei^tcd sum 



determined using the oorrespt ndlng 
of pi evious orthogonal functions, as 



where a is an m X m matrix 
Orthogonal functions from 
a is used thus only the non-: 
used to store a. 



c :)ntaining for Gram-Schmidt 
thp previous orthogonal functions 
zfcro elements are stored thereby 



10 The first Orthogonal functior 



(«)w,(n) 



A time correlation function is defmed as: 
15 and ii can be shown that : 



The weight Omr, relating v^m to Wr can be computed by : 



a«r = — 



TV, 



20 



Since determining the 
arc used to determined the 
detcnnined as follows: 



25 



The weights gm of the onh agonal series is: 



C{n) 
Dim. m) 



where 



ighls used to derive the 
Only half of the matrix 
reducing the memory 



is set equal to tfie first candidate pair (the CDF). 

(8) 



candidate function 
given by: 

(7) 



(9) 
(10) 



(11) 



orthJgonal functions Wn(n) is not re< [uired, as the weights gm 
vcights atn of the non-orthogon i) fimctions» TK^fV) can be 



(12) 



(13) 



12 
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20 



It can be shown that the wcijts relating the candidate functions 
are given by: 



where 



If Ae CDF is less than 1 Vtz 
10 deiermincd in step 2 1 2 tiien 
SCDF is fit as be second 



nd die motor speed is less 
le CDF Is fit as the first 
candidate frequency pair in the 



Afker the first candidate 1ifeq|iency paii(s) 
step 214 or step 216, the 
IS set in step 214 or step 216 
The method 200 attempts to 
functions and the sampted 



imekn 



ad 



The mean sqi»red error is g 



25 



Thctenn gWrnin) is then 
candidate frequency 
the sampled current wave is 
squared error reduction for 
frequency is determined 
reduction for that frequenc; 



an 1 



Each candidate frequency 
if the SCDF is used)hi 



the set 



raO 



i=m+l, m+2, 



Than 2% rated base spejsd as 
candic ate frequency pair and the 
or hogonal set in step 2 1 6. 



has been fit in 
squared error (MSE) betwce 
die sampled current wave is 
minimize the mean squared 
wave. 



currents 



ven by: 



ean squared error reduction 
pair to^e orthogonal set. Since each 
composed of a sine and a 
^h the sine and cosine 
combined to form the 



fin- 



dair in the set is fitted as the 
orthogonal pairs in step 22C 



(H) 

to the sampled signal 
(15) 

(16) 



the orthogonal pairs in either 
1 the orthogonal functions 
Icicrmined in step 218. 
between the orthogonal 



en 3r 



(8) 



caiisedl 



by adding the mth 
frequency component of 
cosiae component, ^e mean 
candid stte of a particular 
overall piean squared error 



second 



frequency pair (or third 
in a manner similar to the 



13 
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fltdngs in steps 216 and 2 14 
(or CDF and SCDF) has alrca|jy 
frequency pairs. The means 
onhogDnal sets resulting firon 



10 



15 



20 



Note that since the frequency 
been fit it will not be in the 
! quared error reduction for 
step 220 is detennined in step|222 



After all of the candidate 
pair if the SCDF was used) ii 
reduction detennined, a 
the highest reduction in the 
Alternatively, all candidate 
of the highest mean squared 
candidate frequency pair at z 
the method 200. 



25 



free uency | 



pairs have been tried as 
the orthogonal set and the rest 
predetermined number 

squared error are selected 
providing a reduction 
jrror reduction may be selected 
time reduces ifae number 



mean! 



pursj 



The selected candidate frequ »icy 
that are CDF harmonics. Tlis 
specific amounts. The remaining 
eh orthogonal pairs in step 2 28 



c srresponding to the CDF 
emaining candidate 
each! additional pair In the 



he second pair (or third 
Itingmcan squared error 
of candidate [Brequency pairs providing 
n step 224. 
within a predetermined range 

Selecting more than one 
of coinputations performed by 



pairs are examined in step 
may be done by looking for 
ft selected candidate 



226 to remove any pairs 
requencies separated by 
freque icy pairs are inserted into 



m( t 



Steps 220 to 228 are repeal d until a specified criteria is 
230. The specified criteria may include fitting a determinejl 
pairs, fitting enough candidi nc pairs to reduce the mean si 
threshold, continue to fit candidate pairs until there is no 
squared error. 



The remaining pairs arc ass imed to be motor speed 
be generated from the rema ning selected candidate 
them to a knovm mathemat cal model relating the induced jfrequcncies 



30 



speed in step 232. An exer iplary mathematical model 
'The Nature of Induction ^ iachines," P. Alger; Gordon an 



Fig- 3 is a system diagram 
includes a controller 304 il 



showing the speed estimator 1 1 
at intercormects a 

14 



as detennined in step 
number of candidate 
error to a certain 
rebuction in the mean 



qiared< 



harmojiics. A speed estimate can 
freque icy pairs by comparing 
to the motor 
suitbble for use is described in 
i Breach, New Work, 1965. 



The speed estimator 1 14 
sampling me ::hanism 306, a connelation 
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mechanism 308. a fitting mec lanism 
controller Interfecc 324, a CE(F 
mechanism 3 1 6, a pair comp^son 
mechanism 320 and a 



The conrroUer 304 coordinate s 
facilitation of data movemen 



The vector controller 
10 indication that the speed 
motor drive coimections 1 
constant CDF for a predetermined period 
is pioduced by the vector 
estimator 1 14 by the vector 



interfa4c 324 interfaces with the vecio 
estimator 1 14 is going to measure 
12|thus allowing the vector 

of time. The CDI 
cobtroller 102 during this time is 
lontroller interfisice 324. 



A sensor interface 300 Is 
convener 302. The sensor 
the motor drive connectiom 
measured current wave 
20 converter 302 for conversicfi 
the sampled current wave 
smaller than the CDF. 



CO mected t 



ma^ 



The A/D converter 302 
25 data. If the current wave 
the current wave need not 



The sampling mechanism 
samples therefrom to 
30 the speed of the motor 1 0^ 



The correlation mechanisip 
samples and determines a 



310, an MSE mechanisfi 
analysis mechanism 3 12, an 
mechanism 318, a 
322 



314, a vector 
4SE reduction 
haimohics recognition 



processing between these i 
and now control. 



mc dules including 



to the sampling 
i iter&ce 300 is in communicatibn 
1 12 to receive the current wav 
at the sensor interface 30C 
from an analog signal to a 
rhe key information In the 



reaived 



be any known A/D converter 
riceived from the sensors 122 i 
K processed by the A/D 



conveiter 



m receives the sampled currei t 
prod ice the sampled current wave t|at 



controller 1 02 to send an 
t he cuirent wave in the 
controller 102 to send a 

(or SCDF and CDF) that 
eceived at die speed 



mecliajdsm 306 through an A/D 
with the sensors 122 at 
5 measured therefrom. The 
is provided to the A/D 
signal, thus forming 
cunlcnt wave is about 60db 



digital; 



::apable of handling such 
is already a digital signal then 
302. 



wave and extracts N 
will be used to estimate 



cunent 



308 receives the sampled 
tonelation between the sampleb 



IS 



wave consisting of the N 
cunent wave and a set of 
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10 



15 



20 



candidate frequency pairs thai are produced by 
current wave may be repicsen t 
frequency pairs. 



The CDF analysis mcchanisn 
controller 102 is less than IH 



thecontr 
d as a weired sum of a 



Uerb04. The sampled 
sub! et of candidat 



3 1 2 determines whether the 
: and the motor speed is less 



indicating low speed operatic i of the motor 104. If the abov^ 
SCDF is used in addition to I le CDF. That is, the CDF anal; 
classifies a previous motor si eed or the speed reference 
nomital. 



th( 



The fitting mechanism 3 JO 
which may also be used to 
fit in the orthogonal set is 
second pair in the orthogona 
weights mechanism 326 and 
weights mechanism 328 
values of the frequency pairi 
v^i^ts for the candidate 



f ts the candidate frequency pair i in an orthogonal set, 
rc present the sampled wave cum nt The first pair that is 
CDF. If the SCDF was used t len the SCDF is fit as *e 
set The fitting mechanism 3 1 0 comprises a frequency 
an orthogonal weigjiis mechan sm 328. The orthogonal 
det4rmines the wcigjits for the orth agonal set based on the 
. The frequency weights mech anism 326 deteimines 
frdquency pairs based on the orthi agonal weights. 



25 



30 



The MSE mechanism 3 1 4 
set and the sampled current 
mean squared error reductic a 
the orthogonal set. The paii 
provide the highest mean 
and the orlhogonal set whei 
comparison mechanism 311 
frequency pairs that provid : 
all candidate pairs providing 
range of the highest value, 



The harmonics recognition 



thm 



from the vector 
2% rated base speed, 
holds true then the 
^sis mechanism 312 
108 (as being eiHier low or 



d ^ermines the mean squared err 
bvave. The MSE reduction 
when additional candidate 
comparison mechanism 318 
sduared error reduction between 
they are added to the oi 
may select a predeteimined 
die highest mean squared errolr 
a mean squared error reductic n 



mechanism 320 determines the se candidate frequency 
pairs from the selected can Hdatc frequency pairs determii ed by the pair comparison 



)r between the orthogonal 
mecjhanisni 316 determines the 
frdquency pairs arc fit into 
d rtermlnes those pairs that 
he sampled current wave 
•ithogor al set The pair 
ntlmber of candidate 
reduction or may select 
within a determined 



16 
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10 



15 



20 



hamw nics of the CDF and Tcmoves i lem from the selected 



mechanism 318 that are I 
candidate frequency pairs. 

The speed estimation mechan sm 322 estimates the speed of 
5 the selected candidate ftequei cy pairs remaining after the 
removed. The speed estimati )n is based on a comparison 
candidate frequency pairs wii \ a known mathematical model 
frequencies in the motor driv : connections 1 12 with the mot >r 



;6012 



Fig. 4 depicts a frequency 
matched for a speed estimate 
the A/D converter 302 is 
this example, the CDF is 
with resulting unwanted odd 
amplitude are not shown 
shown at 25Hz that indicatec 
25Hz/CDF(orCDF 
where the example rated 
are shown to either side 403 
harmonics are shown at 1 20 
that the second frequency, 
consequently, the frequency! may 
determination. 



spdctrum graph of the measured 

by the speed estimator 1 14, 
sho{vn 406 with reference to the 
and drives the motor 104 at 
harmonics 404 (even and 
forptarity). A speed-induced sub 
a rough estimate of motor 
ind SCDF) (60Hz) x rated base 

speed is 1800RPM. Other 
of the sub-harmonic, Exampl* 
■fe and 180Hz 405 separated b] 
ipOHz, coincides with a hannon c 
not be used in the final 



Ibass 



23 



Figs. 5A-B illustrate a speei I 
embodiment of the present 
corresponding with step 20! 
will not be discussed in det ul 



30 



After tihe mean squared 
and SCDF) and the sample I 
speed-induced sub-harmon c 



he motor 104 based on 
harmonics have been 
remaining selected 
relating the induced 
speed. 



CI^F 
of he 



■ spe td 



estimation m&boA SOO accorc ing to a second 
nvcntion. Steps 502 to 5 1 8 of he method 500 
to 21 8 of the method 200 of tl e first embodiment and 



cprrcnt wave, filtered and 
sampling sensitivity of 
vertical (amplitude) axis. In 
rated base speed 401 
otheii harmonics of similar 
larmonic frequency 402 is 
as: 

speed = 750 rpm 
unstated noise frequencies 
of two speed-induced 
2 times the CDF. Note 
of tlie CDF; 
estimate 



spsedi 



err^r between the orthogonal set 
current wave has been 
frequencies is performed 



c( mtaining the CDF (or CDF 
deterai ined, a search for possible 
betjveen the CDF (or SCDF) 



17 
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and OM2 in step 520. 
induced firequencies 



10 



15 



20 



Sub-harinonic frequencies arc those fio juencies 
hannonlc 5 but have a lower frequency ttlan 



25 



Each speed-induced sub-hanr Dnic determined in step 520 is 
relating the induced frcsqucnc bs with the estimated speed, 
comparison is to locate two s< parate frequency regions when : 
induced harmonics would ha) e to be in order to make a pattf m 
model in step 522. These frei lucncy regions correspond wit^ 
correlation in step 510. Thes t regions may be equivalent to 
FFT bin lower than the Identi Bed sub-harmonics. 



compared to a model 
purpose of this 
corresponding speed 

match wich the speed 
the FFT bins used m the 
FFT bin higher and I 



Any frequency harmonics 
are then located in step 524 
be used to filter and/or 
to avoid the lime wasted 
of Ae induced harmonics). 



of die < 



sub-hannonics in each i 
During a search of each 
discri ninate against other noise cairibrs 
deu cting and processing erroneous 



The speed-induced frequenc|r 
harmonic frequencies 
two tnncs the CDF away in 
by the smallest harmonic 



identified! 



frequency region 
region[variou$ techniques may 
using learned methods 
signals (e.g. an inversion 



frequency 



harmonics found are examinep 
corresbond to a second harmonic 

tep 526 or two harmonic ft«qu 
contto) signal (e.g. CDF or 



of the < 



hi step 528 it is determined jvhcther there were any frequcti zy 
526, If there was ai least on 5 frequency pair was dctermin^ 
that has the greatest separat on fit>m the location of (mown 
harmonics) is selected hi st^ 530 to avoid unwanted motoi 

st^p 



frequencies. The selected fi equency pair is compared in 
harmonic speed model to es timaic ihc motor speed. If therj^ 
speed estimate, the estimate 5 are averaged to result in a 



30 conm>ller 104 in step 534. 

If there were no frequency 
536 if other sub-hannonics 



fmd 



>airs determined in step 528 
were found. If other 

IS 



that are speed- 
thcCDF(orSCDF). 



to determine if the 

whose frequency is 
^ncies that are separated 
3DF+ SCDF). 



pairs determined in step 
then the fiiequency pair 
larmonics (c-g. CDF 
rather dian speed related 
532 to the known 
is more than one valid 
value to drive the vector 



ih^n it is determined m step 
sub-harmimics were found then steps 



04-F8lrl3 09:4Sain From-KIRBY EADES GALE BAKER 



+16132370045 



T-704 P. 032/065 F-051 



522 to 528 are repeated. If no 
estimate is produced in step 5 



Matching the two fteqtiencies 
s speed-induced harmonics 
with known hannonics may 
sophisticated pattern matchin 



froii 



b5 



Fig 6 is a system diagram 
ID embodiment The speed 
sampling mechanism 606, a 
MSE mechanism 614, a 
612, a sub-harmonics 
hannonics identification 
15 and a speed estimation 



sh<|wing the speed estimator 114 
114 includes a 
i omlation mechanism 608, a 
controller interface 624, a 
mechaldsm 615, a region detcnninaticfi 
mec hanism 610 a correspondhig 
mech mism 622. 



estii aator 



vect 



The conn-oiler 604 coordinai ss 
facilitation of data movemer t 



20 The vectOT conttoDer 
indication tiiat die speed 
motor drive connections 1 1 
constant CDF for a 
is produced by the vector 

25 estimator 1 14 by the vector 



A sensor inter&ce 600 is 
convener 602. The sensor 
the motor drive connection! 
30 measured current wave 

converter 602 for conversi({n 
the sampled current wave, 
smaller than the CDF. 



other sul>-hannonics were i 
8. 



fow d then no new speed 



means! 



as in step 526 is not the only 
unwanted harmonics, 
discarded. These steps may 
algorithms and a speed trackukg 



of separating 
Frequehcy pairs that coincide 
replaced whh more 
routine. 



tei 



a xording to tfie second 
controller S04 that interconnects a 
fr ting mechanism 610, an 
analysis mechanism 
mechanism 618, a 
luencies mechanism 626 



processing between diese monies including 
and flow control. 



interfej^e 624 interfaces with the vectcjr 
es^ator 1 14 is going to measure 
thus allowing the vector 
predeterhined period of time. The CD ' 
cc niroller 102 during this time 
controller interface 624, 



ccnnectedl 



; mccha lism 



to the sampling 
nter&ce 600 is in communicatil^n 
1 12 to receive the current 
reciivcd at the sensor interface 30( 
from an analog signal to a di; 
The key information in the curfent 



controller 102 to send an 
he current wave in the 
controller 102 to send a 

(or SCDF and CDF) that 
isteceived at the speed 



di; ;ital 



606 through an A/D 
with the sensors 122 at 
: measured therefrom. The 
is provided to the AID 

signal, thus forming 
wave is about 60db 



19 
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The A/D converter 602 may b 
data. If the cuirem wave rece 
the current wave need not be 



5 any known A/D converter 
ved from the sensors 122 is 
iroccssed by The A/D converter 



The sampling mechanism 60i 
samples therefrom to produc^ the sampled 
the speed of the motor 1 04. 



6)8 1 



10 The correlation mechanism 
sampks and determines a 
candidate frequency pairs 
current wave may be 
frequency pairs. 

15 

The CDF analysis mechanisrji 
controller 102 Is less than 
indicating low speed operation 
SCDF is used in addition to 
20 classifies a previous motor 
normal. 



receives the sampled 
coi^clation between the sampled 
are produced by the controllei 
represented as a weighted sum of a 



The fitting mechanism 610 
which may also be used to 

25 fit in the orthogonal set is th 
second pair in the onhogons i 
weights mechanism 630 anc 
weights mechanism 628 
values of the frequency pair^ 

30 weights for the candidate 



The MSE mechanism 614 
set and the sampled current 



cadable of handling such 
air iady a digital signal ihcn 
6302. 



receives the sampled current vf avc 
current wave that 



currenll wave consisting of the N 
cprrent wave and a set of 
604. The sampled 
of candidate 



sulset 



CDF 



dan 



612 determines whether the 
and the motor speed is less 
of the motor 104. If the abo^ 
he CDF» That is, the CDF analysis 
speed or the speed reference 108 



litsi 



' pan 3 



the candidate fi^uency 
present the sampled wave 
CDF. If the SCDF was used 
set. The fitting mechanism 6 
an orthogonal weights mechaif sm 
determines the weights for the 

J. The frequency weights mecl 
fin quency pairs based on the orth >gonal weights. 



(ttermines the mean squared ertpr between the orthogonal 
wave. 



and extracts N 
wiW be used to estimate 



from the vector 
2% rated base speed, 
e holds true then the 
is mechanism 312 
as being either low or 



in an orthogonal set, 
cuiT^m. The first pair that is 
1 hen the SCDF is fit as the 
0 comprises a frequency 
628. The (vthogonal 
orttlogonal set based on the 
anisro 630 determines 



20 
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10 



The sub-hannonics mechanisi^ 
current wave between OHz 



616 locations all specd-induc ;d sub-hannonics in the 
anji the CDF (or SCDF and CDF) 



5 The region determination me( 
harmonics with a harmonic s 



hanistn 618 compares one oft) e determined sub- 
eed model to identilV frequenc ^ search regions. 



The harmonics identification 
sub-haimonics in the identifiild 



ncchanism 620 locates speed- 
frequency search regions in 



The corresponding fkequencii s 
regions determined by the re|ion 
corresponding frequency in 
CDF between the two 
15 the smallest harmonic of the 
one frequency pair is 
from other harmonics is 



mechanism 626 

detenninaticm mechanism 
second region where there ii 
1, or two harmonic frequi 
jontrol signal (e.g. CDF or CD 
frequency pair having 
pro^bded to the speed estimation 



tie! 



I freque icies, 



located then the j 



The speed estimation 
20 the selected candidate 

removed The speed 

candidate frequency pairs 

frequencies in the motor 

corresponding frequency pairs 
25 from each of diese pairs is aferaged. 



with 
dri>e 



lan{ uagc. 



Embodiments of the 
computer programming 
in a procedural programmin 
30 "C++")' Further cmbodimejits 
programmed hardware 
hardware and software components. 



iduced harmonics of the 
the current wave. 



determines frequencies in one of the 
18 that have a 
is a spread of two times the 
lenties that are separated by 
?+SCDF. If more than 
the greatest separation 
m< chanism 622. 



of the 



CDF 



mechaiiism 622 estimates the speed 
ftequdncy pairs remaining after the 
estimaqion is based on a comparison o: 
a known mathematical mode 1 
connections 1 12 wi0i the mo 
have been located then the 



present invention may be implementec 
For example, 
I language (e.g. *'C") or an 
of the invention may be 
clenri|snts, other related componentsJor 



motor 104 based on 
harmonics have been 
the remaining selected 
relating the induced 
or speed. Ifnniliiple 
estimate obtamed 



sjieed 



in any conventional 
embodimehts may be implemented 
objcjct oriented language (e.g. 
im[|[ementjed as pre- 

as a combination of 



21 



I impleme ited 



; embod es all 



art 



Embodiments can be 
computer system. Suchii 
fixed either on a tangible 
diskette, CD-ROM, ROM, or 
5 modem or other interface 
network over a medium 
or electrical communications 
techniques (e.g. microwave 
computer instructions 
10 herein. Those skilled in the 
be written in a number 02 
architectures or 
any memory device, such as 
devices, and may be transmit 
15 optical, infrared, microwave, 
such a computer program 
accompanying printed or 
preloaded with a computer s; 
from a server over the 
20 embodiments of the inventic^ 
(e.g. a computer program 
embodiments of the invenii<|n 
software (eg. a 



• operating sys terns. 



L computer pi sgram 



25 It is apparent to one skilled 
from the specific embodimcbts 
the spirit and scope of the invention 



pre duct 



as a computer program 
implehentation may include a series 
J med um, such as a computer rcadab 
Sxed disk) or transmittable to a 
J devi ce, such as a communications 
Thelmedium may be either a tangible 
iocs) or a medium implemente i 
ijifr^red or other transmission te 
or part of the fimctionalit i 
should appreciate that such c )mputer 
if pro|;Rnnming languages for use with 

Furthermore, such instrufctions 
©miconductor, magnetic, opti4' 
Ded tising any communications 
or other transmission tecfanolo, ^es 
product may be distributed as a re 
cle ^onic documentation (e.g. shrink 
astern (e.g., on system ROM 
: (e.g., the Internet or World 
may be implemented as a 
pil>duct) and hardware (termed m 
may be implemented as entir 
product).** 



n the art that numerous 
described herein may be 



for use with a 
3f computer instructions 
e medium (e.g. a 
c mputer system, via a 
^apter cormected to a 
medium (e.g, optical 
with wireless 
:hnlques). The series of 
previously described 

instructions can 
many computer 

may be stored in 
or other memory 
echnology, such as 

It is expected that 
novable medium with 
wrapped software), 
disk)> or distributed 
>/idcWeb). Some 

of both software 
chanisms). Still other 
ily hardware, or entirely 



or fixed < 



cor ibination ( 



modifidations and departures 
mjde wtthout departing from 



22 



